Pauli Spin Blockade in Carbon Nanotube Double Quantum Dots 
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We report Pauli spin blockade in an impurity denned carbon nanotube double quantum dot. We 
observe a pronounced current suppression for negative source-drain bias voltages which is investi- 
gated for both symmetric and asymmetric coupling of the quantum dots to the leads. The measured 
differential conductance agrees well with a theoretical model of a double quantum dot system in the 
spin-blockade regime which allows us to estimate the occupation probabilities of the relevant singlet 
and triplet states. This work shows that effective spin-to-charge conversion in nanotube quantum 
dots is feasible and opens the possibility of single-spin readout in a material that is not limited by 
hyperfine interaction with nuclear spins. 

PACS numbers: 73.23.-b, 73.63.-b, 75.10.Jm, 85.35.-p 



I. INTRODUCTION 

The electron spin is a natural two-level system and 
therefore attractive as a quantum bit in quantum infor- 
mation processing schemes. Spin qubits defined in quan- 
tum dots are of particular interest because of the possi- 
bility to isolate, manipulate and measure single spins 
Much of the attraction of spin qubits in quantum dots 
is also related to the relatively long time over which a 
superposition of opposite spin states of a single electron 
remains coherent. This long spin coherence time is a di- 
rect result of the electron's small magnetic moment which 
ensures it couples only weakly to its environment. The 
inevitable problem this poses for the read-out of a single 
spin is elegantly overcome in double quantum dot sys- 
tems by converting the spin information to a charge state 
using the phenomenon of Pauli spin blockade which oc- 
curs when certain transitions between two quantum dots 
are forbidden by spin selection rules 0, Q as illustrated 
in Fig. 1. Significant experimental effort on spin qubits 
defined in quantum dots in GaAs/AlGaAs heterostruc- 
tures has already resulted in the demonstration of driven 
coherent oscillations of single electron spins [J] and the 
coherent exchange of two electron spins in a square-root- 
of-swap quantum operation 

While significant as a proof-of-principle that single spin 
manipulation and read-out in a solid-state environment 
is feasible, these experiments also demonstrated that the 
spin coherence time in these devices is limited by hyper- 
fine interactions with the Ga and As nuclei 0, H, H 0| • 
As a result, the number of coherent single-spin rotations 
that can presently be observed within the spin coherence 
time in GaAs based quantum dots is several orders of 
magnitude below the typical figure-of-merit of 10 4 quan- 
tum operations for a fault tolerant quantum computer. 

There is therefore a strong incentive to develop spin 
qubits in materials in which hyperfine interactions are 
much reduced or absent altogether. Carbon based mate- 



rials such as carbon nanotubes or graphene are excellent 
candidates in this respect. Due to the absence of hy- 
perfine coupling in the dominant 12 C isotope, the spin 
coherence times are expected to be very long 8, 9] while 
the recent observation of spin-orbit interaction in nan- 
otube quantum dots suggests the possibility of electrical 

control of the spin states 0, m d, m. in this work, 

we show that spin blockade is readily observed in weakly 
coupled carbon nanotube quantum dots even for many 
electrons on the nanotube and for temperatures of order 
one Kelvin. We therefore conclude that effective spin- 
to-charge conversion in carbon nanotube quantum dots 
is feasible and that single-spin manipulation and readout 
in nanotubes is a promising and realistic prospect. 



II. NANOTUBE DEVICES 

The device we consider consists of an individual carbon 
nanotube filled with Sc@Cg2 molecules (that is a Sc atom 
inside a Cs2 cage) contacted by electron-beam defined 
palladium source (S) and drain (D) electrodes that are 
separated by 300 nm. The degenerately doped Si/SiC>2 
substrate (300 nm oxide) is used as back gate. The study 
of carbon nanotubes filled with Sc@Cs2 is motivated by 
the long spin coherence times of the unpaired spins on 
the encapsulated Sc atoms as observed by us in ensemble 
measurements [l3j . 

Measurements on twelve different devices suggest that 
even though the presence of Sc@C§2 may lead to ob- 
servable bandstructure modification or charge transfer 
(doping) between the nanotubes and Sc@Cs2, the low- 
temperature transport properties are remarkably similar 
to those of empty nanotubes [3] . These findings are con- 
sistent with recent transport experiments on nanotubes 
with encapsulated Ceo molecules as well as with 

recent density functional theory calculations of spin in- 
teractions of chains of Sc@Cs2 inside carbon nanotubes 
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FIG. 1: (a) Schematics of the electrochemical potentials of the relevant one and two-electron states of a double quantum dot 
in the absence of a tunnel coupling between them. The level offset between the single-particle states on both dots is given by 
Se. The on-site charging energy and electrostatic coupling energy are shown as U and U' , respectively, (b) When the tunnel 
coupling t between the quantum dots is significant, the S(l,l) and S(0,2) singlet states hybridize to form molecular bonding 
(Si) and anti-bonding (Sa) singlet states which separate from the T(l,l) triplet states, (c) When a negative bias voltage is 
applied and one of the T(l,l) triplet states becomes occupied, the (1, 1) — > (0,2) transition is not allowed because, by virtue 
of the Pauli principle, the (0,2) state has to be a spin singlet and any further current flow is blocked, (d) For opposite bias 
conditions the (0, 2) — » (1, 1) transition through the singlet states is allowed. Note that for finite detuning, the Si and S2 singlet 
states (as well as the triplet states) are still extended molecular states but are dominated by the (0,2) and (1,1) charge states 
in the way indicated by the schematics. 



|17| . The conclusions of our present work will therefore 
apply equally well to empty carbon nanotubes. 



III. RESULTS AND DISCUSSION 

The linear-response conductance of the device as a 
function of gate voltage measured at temperature T = 1.4 
K is shown in Fig. 2(b) which displays a series of irregular 
conductance peaks. A variation in peak height and spac- 
ing is common for carbon nanotube quantum dots and 
is generally attributed to scattering by defects along the 
nanotube. As demonstrated by a combination of scan- 
ning probe and transport experiments [18], structural de- 
fects in particular lead to resonant electron scattering in 
which the defects can be transparent or opaque depend- 
ing on the electrochemical potential (and hence the gate 
voltage). Defects could be introduced in our nanotubes 
by, for example, the acid treatment used to purify the 
samples from magnetic impurities. However, defects are 
also commonly observed in as-grown nanotubes such as 
those made by chemical vapor deposition and for which 
typical scattering lengths of ~ 150 nm were reported p^j . 
Low-temperature transport experiments on a number of 
our nanotube devices of the geometry studied here indi- 
cate that most consist of a series of (two or three) quan- 
tum dots. 

Double quantum dots in carbon nanotubes in which 
the tunnel barrier between the dots is due to a defect in 
the nanotube are relatively common and have been stud- 
ied previously in, e.g., Refs. [H, [2(|. Typically, the dou- 
ble quantum dot is identified by the characteristic hon- 
eycomb pattern that develops as two independent gate 
electrodes (each coupled to a different quantum dot) are 
varied [2l| . In our device we only have the ability to vary 
a single gate electrode (the back gate). Nevertheless, a 
detailed understanding of the physical phenomena under- 
lying the electronic transport properties can be obtained 



by studying the differential conductance (dl/dV) as a 
function of gate (V g ) and bias voltage (V s d), as demon- 
strated for double quantum dots by Ono et al Q. Here 
we focus on two pairs of peaks around V g = —1.95 V and 
V g = —3.65 V of which the the differential conductance 
is shown in Figs. 2(c) and 2(d), respectively. The most 
striking feature of the dl/dV plots is the pronounced 
asymmetry in the bias voltage and the current suppres- 
sion and appearance of negative differential conductance 
for negative V s d (dark blue regions). Also note that while 
Fig. 2(c) is approximately symmetric in gate voltage, this 
symmetry is clearly broken in Fig. 2(d). 

The lack of periodicity in the linear-response conduc- 
tance, the pronounced negative differential conductance, 
as well as the striking difference in the slopes of adja- 
cent Coulomb diamonds (most apparent in Fig. 2(d)) are 
clearly at odds with a model of a single quantum dot. As 
we will show below, these features can be explained well 
by invoking Pauli spin blockade in a double quantum dot. 
The simplest model of two coupled single-level quantum 
dots as introduced in Fig. 1 suffices to explain our mea- 
surements. This might seem surprising given that the 
nanotubes will contain many electrons but is justified 
by the large single-electron level spacing AE (see below) 
and the simple even-odd shell filling of carbon nanotube 
quantum dots [22j • Recent work on GaAs double dot 
systems containing more than two electrons (up to ~ 10) 
could also be explained in terms of effective single-level 
quantum dots [23|, H3] ■ 



A. Double quantum dot Hamiltonian 

To quantitatively compare the measurements in Fig. 2 
with an interpretation in terms of spin blockade, we have 
used a many-body density matrix approach to calculate 
the current and population numbers of the eigenstates 
of two coupled single-level quantum dots. The double 
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FIG. 2: (a) Schematic level diagrams for different transport 
regimes corresponding to the measurements of panel (c), see 
symbols, (b) Linear-response conductance of the device mea- 
sured at T = 1.4 K. (c) Color-scale plot of the differential 
conductance (dl/dV) as a function of source-drain bias volt- 
age (V s d) and gate voltage (V g ) for the gate region indicated 
by the arrows in panel (b). Dark blue corresponds to neg- 
ative dl/dV. The ordered pairs (n,m) indicate the effective 
electron occupancy in each quantum dot. (d) Differential con- 
ductance for a different region of V g , showing a pronounced 
asymmetry in both bias and gate voltage. 

quantum dot (DQD) and the leads are modelled by the 
Hamiltonian H = Hl +Hr+Hdqd +T~Lt, where Ti-L(R) 
models the left (right) lead in a free electron like approx- 
imation. The double quantum dot is modelled by 



Hdqd = f e ncJ d) na d, ia + U n n ^n„i ) 

n=A,B cr ' 

+ U'(nA-\ + n Ai )(n B ] + n Bi ) 

+ J2( td AjB. + H.c). (1) 

cr 

In the model we have two levels, £a<j and £_b ct , where 
cr =t, I is the the electron spin. Here, d' A , g s c (dA(s)o-) 
creates (annihilates) an electron in quantum dotA(B) 
with spin a. The on-site Coulomb charging energy and 
electrostatic coupling energy are denoted by U and U' , 
respectively, whereas t is the tunneling rate between the 
quantum dots, see also Fig. 1. The last term in TL 
accounts for the tunneling between the leads and the 
DQD. This model neglects spin-orbit and hyperfme in- 
teractions. 

The Hamiltonian "Hdqd is transformed into diagonal 
form, e.g. H DQD = Y / Nn E Nn\N,n)(N,n\, where E Nn 
is the energy for the eigenstate \N,n) with N electrons, 
where n is a state label (in our model there are 16 eigen- 



states, such that N = 0, n = 1, N = 1, n = 1, . . . ,4, 
N = 2, n = 1,...,6, N = 3, n = 1,...,4, and 
N = 4, n = 1). This enables a many-body density 
matrix approach for calculating the population number 
probabilities Pjy n of the corresponding state \N,n), to 
the first order approximation with respect to the cou- 
pling r L ( R ) to the left (right) lead. This order of ap- 
proximation is based on that only the diagonal transi- 
tions \N,n) — > \N, n) are included, while effects from 
off-diagonal transitions such as \N,n) — > \N,n') would 
require a higher order expansion of the rate equations 
for the population number probabilities. Our calcula- 
tion provides the dynamics of the populations numbers 
as function of the bias voltage and equilibrium electro- 
chemical potential. Knowledge of the population number 
probabilities enables calculation of the current and dif- 
ferential conductance through the system, using standard 
techniques. The approach provides complete knowledge 
of the matrix elements (N, ii\Ht\N±1, to), for transitions 
between states differing in electron number by one. This 
information allows for detailed analysis of which states 
are involved in the conductance and which state(s) are 
responsible for spin and Coulomb blockade. The method 
is more thoroughly described in Ref. 



B. Double quantum dot model parameters 

We start by comparing the model predictions with the 
measurements of Fig. 2(c). In the model, and in the anal- 
ysis below, we assume that the charging energy U and the 
capacitive coupling to the gate electrode is the same for 
both quantum dots. This approximation is justified by 
the symmetry (in V g ) of the data in Fig. 2(c). The re- 
sult of the model calculations, using an appropriate set 
of parameters which we will discuss below, is shown in 
Fig. 3(a). The calculated differential conductance is in 
good agreement with the experimental data of Fig. 2(c) 
for the low bias regime while differences are observed at 
higher bias voltages. Figure 4 illustrates how the differen- 
tial conductance (evaluated in Fig. 4 for the electrochem- 
ical potential corresponding to the dotted vertical line in 
Fig. 3(a)) is related to the occupation probabilities of the 
various one and two-electron states. These calculations 
also show, see Fig. 4(b), that the observed current sup- 
pression for negative V s d is indeed the result of a near 
unity occupation probability of the T(l,l) triplet states. 

The first model parameter to consider is the voltage 
drop AVqd over the tunnel barrier between the two 
quantum dots in the presence of a source-drain bias volt- 
age. The relative voltage drop can be obtained from the 
slopes of the diamonds in the experimental data, marked 
by L and R in Fi g. 2 (c) and from an electrostatic model 
of the device and yields AV QD /V sd - 0.15. As 

a result, the bias voltage V s d acts as a knob that effec- 
tively controls the level offset 5e (or 'detuning') between 
the quantum dots plj . The magnitude of Se can be ex- 
tracted from the excitation line indicated by the triangle 
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FIG. 3: (a) Color-scale plot of the calculated differential con- 
ductance as a function of source-drain bias voltage and gate 
voltage. The model parameters used are &e = 5 meV, U' — 12 
meV, U = 17 meV (both quantum dots) and t — 0.7 meV. 
The relative voltage drop over the tunnel barrier between the 
two quantum dots is AVQD/Vsd ~ 0.15. The temperature is 
set to T = 4 K. The symbols correspond to the schematics 
in Fig. 2. (b) Same as panel (a) for asymmetric coupling of 
the double quantum dot to the leads as described in the main 
text. 



in Figs. 2(c) and 3(a). The position of the triangle in the 
dl/dV plots corresponds to the situation in which (ne- 
glecting the tunnel coupling) the electrochemical poten- 
tials on the left and right quantum dot are aligned with 
the electrochemical potentials of the source and drains 
electrodes, respectively, see right schematic in Fig. 2(a). 
Taking into account the voltage drop between the quan- 
tum dots, this yields 8e ~ 5 meV. When the bias volt- 
age is increased beyond this point, the model calculation 
shows, see Fig. 4(b), that the double dot gets trapped in 
a (1,0) charge state (this excitation line is indicated by 
the lower arrow in Fig. 3(a)). As a result, the current 
decreases again with increasing V s d- In the experiment, 
however, this is not observed. This difference can be un- 
derstood considering that the model does not account for 
inelastic scattering processes. This approach is justified 
in the spin blockade regime since relaxation due to e.g. 
electron-phonon interaction is strongly suppressed if the 
transition involves a spin-flip [26| . On the other hand, 
the (1,0) — > (0,1) transition does not require a spin-flip 
and energy relaxation by phonon emission will be effec- 
tive in the experiment. As a result, an increase in current 
is expected (and indeed observed) along this excitation 
line, as indicated by the arrow in Fig. 2(c). 

A similar reasoning would explain the absence of nega- 
tive dl/dV that is observed in the model for positive V s d- 
For example, the excitation line observed in the experi- 
ment (Fig. 2(c), square, and upper arrow in Fig. 3(a)) 
could be due to a transition between the S2 and Si 
states which does not require a spin-flip either. How- 
ever, since the observed excitation line is exactly parallel 
to the ground state transition, we tentatively attribute 
it to a single-dot excitation, which yields AE ~ 4 meV. 
If we assume the conventional nanotube dispersion rela- 
tion, the level spacing is related to the nanotube length L 
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FIG. 4: (a) The current and dl/dV can be understood from 
the evolution of the two singlet states (Si, red and S2, ma- 
genta) and the triplet states (green) as a function of V B d- 
The model parameters used here correspond to the situation 
of Fig. 3(a). The diagonal lines indicate the opening of the 
bias window as V a d is increased. The electrochemical poten- 
tial /i = 6.1 meV corresponds to the dotted vertical line in 
Fig. 3(a). (b) Calculated occupation probabilities for the var- 
ious one and two-electron states as a function of V s d- (c) 
Peaks in the current and dl/dV are observed when transi- 
tions between the one-electron and two-electron states are 
energetically accessible. 



through AE = Uvf/^L, where h is Planck's constant and 
v F = 8.1 x 10 5 m/s is the Fermi velocity [27j]. For AE ~ 4 
meV this yields L = 200 nm which, given the source- 
drain separation of 300 nm, would imply two quantum 
dots of similar size. 

The electrostatic coupling energy U' can be obtained 
from the size of the main diamond (positive half) in 
Fig. 2(c) from which we obtain U' ~ 10 — 15 meV [2l| . 
The charging energy on the individual quantum dots is 
more difficult to determine exactly. The most satisfac- 
tory correspondence between the data and model is ob- 
tained for U ~ 17 meV for both quantum dots, consistent 
with dot lengths of order ~ 100 - 200 nm [22]. We veri- 
fied that the conclusions of our work are not sensitive to 
the exact value of U. 

The ratio of the tunnel coupling to the level offset is 
crucial for the observation of spin blockade and t/Se <C 1 
must be satisfied for it to be clearly observed To 
appreciate this, consider that a current flow in the double 
dot involves transitions between the one-electron state 
and the two-electron singlet or triplet states. The one- 
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FIG. 5: (a) Measured current as a function of V s d at gate 
voltage V g — —2.02 V, corresponding to the dotted vertical 
line in Fig. 2(c). The red oval indicated the regime of Pauli 
spin blockade with measured leakage current of ~ 100 pA. 
(b) Calculated current traces as a function of electrochemical 
potential at V s d = —5.4 mV for the model parameters of 
Fig. 3(a). The tunnel couplings are t = 0.7 meV (black) and 
t = 0.3 meV (red). Inset: maximum occupation probability 
of the triplet states. 

electron state is a superposition of the (0,1) and (1,0) 
charge states i.e. of the kind a|0, 1) + /3|1,0) and the 
probability of finding the electron in the energetically 
excited (1,0) state depends directly on t/Se [2l|. Because 
in the spin blockade regime, transferring an electron to 
the drain electrode involves a transition from a triplet 
state to the (1,0) charge state, a current flow will be 
strongly suppressed if ft -C a. Note that this is not the 
case for a current mediated by the singlet states which 
follows the sequence (l,l)-> (0,2)-> (0,l)-> (1,1). 

In our experiments, the tunnel coupling can be ex- 
tracted from the small leakage current of ~ 100 pA that 
is observed in the spin blockade regime, see Fig. 5(a). As 
compared to the ~ 5 nA measured at positive V s d this 
implies a current suppression of a factor of ~ 50. Given 
the previous estimate for Se of approximately 5 meV, the 
model requires t < 0.7 meV to provide a similar suppres- 
sion factor. The observed leakage current also directly 
puts a lower bound of e/I ~ 2 ns on the spin relaxation 
time T\ in carbon nanotube quantum dots. Note however 
that this is likely to be a strong underestimate of the in- 
trinsic spin-flip relaxation time in carbon nanotubes and 
that the observed leakage current can be fully accounted 
for by transitions mediated by the remaining finite occu- 
pation probabilities of the one-electron and two-electron 
singlet states. The strong dependence of the leakage cur- 
rent on the tunnel coupling is illustrated in Fig. 5(b). 

C. Asymmetric coupling to the leads 

While the measurements of Fig. 3(a) are approximately 
symmetric in gate voltage, this is not the case for the dif- 
ferential conductance shown in Fig. 3(b). As compared 
to the symmetric situation, there is a pronounced tilt in 
the slopes of the Coulomb diamonds and strong negative 



differential conductance is only observed along the right 
edge of the Coulomb diamond. In the model these fea- 
tures are reproduced by introducing an asymmetry in the 
potential drop at the source and drain electrodes and in 
the tunnel couplings T s /Td of the leads to the quantum 
dots while keeping other parameters such as U, U' and Se 
identical to the symmetric situation. The tunnel coupling 
is set to t = 1.2 meV. The result for AV s /AV d = 3/2 and 
Ts/r^ = 1/5 is shown in Fig. 3(b) which corresponds well 
with the experimental data, see Fig. 2(d). 

The dependence of the current in the Pauli spin block- 
ade regime on the asymmetry in the tunnel couplings for 
the model used here is described in detail in Ref. [3| . An 
intuitive way to understand the effect of the asymmetry 
on the triplet occupation probability (and therefore on 
the leakage current) is the following: in the spin blockade 
regime as considered above, an electron has a high prob- 
ability to enter the double dot from the source electrode 
to form a T(l,l) triplet state but a low probability to exit 
to the drain (hence the high occupation probability) . An 
asymmetric coupling T s <C T^, on the other hand, has 
the precise opposite effect. The result is a reduction in 
the triplet occupation probability such that negative dif- 
ferential conductance is not observed along the edge of 
the Coulomb diamond in the left part of Fig. 3(b). 

The strong negative differential conductance along the 
Coulomb diamond edge on the right side of Fig. 3(b) 
can be understood considering the (broken) electron-hole 
symmetry of the double dot. Whereas the linear-response 
conductance at point E in Fig. 3(b) corresponds to an 
electron moving from the source to drain electrode, the 
conductance at point H corresponds to a hole moving 
in the opposite direction [2l|. While the asymmetry in 
the tunnel rates (partially) lifts the spin blockade for the 
electron cycle, it enhances the blockade for the hole cycle 
(as it moves in the opposite direction). The result is a 
strong T(l,l) triplet occupation probability and negative 
conductance in the dl/dV. 

Note that several of the excitation lines in the model 
(indicated by the arrows) are also observed in the data. 
The difference in the polarity of the excitation (nega- 
tive versus positive dl/dV in the model and experiment, 
respectively) is attributed to energy relaxation by e.g. 
phonon emission which is not accounted for in the model. 
As discussed above, inelastic scattering will be effective 
in the experiment for transitions that do not require a 
spin- flip. 



IV. CONCLUSIONS AND OUTLOOK 

In conclusion, we present measurements of the differen- 
tial conductance of an impurity-defined carbon nanotube 
double quantum dots showing Pauli spin blockade. The 
measurements arc well described by a theoretical model 
of the device which allows us the estimate the relevant 
singlet and triplet occupation probabilities of the dou- 
ble quantum dot. Since the phenomena of spin blockade 
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enables spin-to-charge conversion in quantum dots, our 
findings present an important step towards single-spin 
read-out and spin qubit operations in carbon-based de- 
vices that are not limited by hyperfine interactions. An 
additional advantage of nanotubes in this respect are the 
large energy scales observed here which compare favor- 
ably to other systems such as lateral GaAs or Si double 
quantum dots. In electron-spin resonance (ESR) experi- 
ments [4[ , which are limited by photon assisted tunneling, 
this would allow considerably larger oscillating fields and 
hence much faster single-spin rotations. In fact, the sig- 
nificantly larger g-factor of carbon nanotubes (g « 2) as 
compared to GaAs devices (g k —0.44) would already 
provide a fivefold gain. 

For future experiments, control of the individual quan- 
tum dots and the tunnel coupling is imperative. Fully 
tunable carbon nanotube quantum dots have already 
been studied by a number of groups but spin blockade 
had not been previously observed [2(| HI) HI, H3|- 
This seems surprising given the clear signatures seen here 
at a relatively high temperature of 1.4 K. This difference 
is likely to be related to the condition t/Se <C 1 which 



must be satisfied for the observation of spin blockade in 
double quantum dots. As our device is rather small (dot 
sizes of order 100-200 nm) and does not make use of 
metal top gates that would add to the overall device ca- 
pacitance, most quantum dot energy scales are nearl y an 
order of magnitude larger than in Refs. [l|| H(| [28ll29|. 
We expect that when the device dimensions and/or the 
tunnel couplings are sufficiently reduced, spin blockade 
will also be observed in top gated carbon nanotube quan- 
tum dots. 
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